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ABSTRACT 



A divergent geostrophic prediction model proposed by Duthie is 
investigated to determine its usefulness in producing realistic fore- 
casts and to determine the feasibility of further research using the 
gradient wind. 

24-hour forecasts were made using the prediction equation, 



at 300 mb, 500 mb, and 850 mb for summer and winter situations, with and 
without term B and with a tuning constant on term D. 

Tables presenting the error fields for these tests are included, 
along with the plotted forecast and verifying maps. These tables and 
illustrations will show that the best results at 500 mb and 300 mb are 
obtained with the B term out and a tuning constant of minus four and 
minus one respectively on the D term. At 850 mb different tuning con- 
stants produced no significant changes in the prognostic charts* 

The results of this investigation show that usable prognostic maps 
are produced, that the model tested has operational potentiality, and 
that further testing is warranted. 
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1 . 



Introduction. 



With the advent of the high speed digital computer, forecasting 
shifted from a highly subjective art accomplished by individuals using 
hand-drawn analyses and various prognostic methods, many of which had no 
theoretical foundation, to a more sophisticated method of solving the 
complicated differential equations governing the motions of the atmos- 
phere. Many theoretical models have since been tested and the proficien- 
cy of the forecaster seems to have improved. However many of these 
models have experienced difficulty in moving troughs and systems with 
sufficient speed. 

The purpose of this paper is to introduce a divergent geostrophic 
model proposed by Duthie [2]. This model allows poleward moving parcels 
to increase their relative vorticity and equatorward moving parcels to 
decrease their relative vorticity, thus allowing for development and fill 
ing. It will also move waves with a speed of the basic zonal current. 

The objectives are to show that this model, which uses the diver- 
gence of the geostrophic wind, will: 

a. numerically converge and produce a recognizable weather 

map; 

b. move waves at a speed consistent with the verifying analyses 
and c. produce results sufficiently accurate to warrant further 

investigation. 

The author wishes to express his appreciation to Professor W. D. 
Duthie of the U. S. Naval Postgraduate School, on whose work this paper 
is based, for his guidance and encouragement in this investigation. In 
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addition, appreciation is also expressed to the personnel of the U. S. 
Fleet Numerical Weather Facility for their assistance in this project. 
Special thanks are extended to Mr, Leo Clarke for assistance in program- 
ming and many helpful suggestions. 
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2. Background. 



The isobaric vorticity equation for frictionless flow is [3] 



* - ^ V 1 V + ~T> 



( 2 . 1 ) 



It can be shown [1] that the divergence of the geostrophic wind 



is 



^7.y^ = - Vy Any f . 



(2.2) 



Assuming 



CO £ = T 

O qP 



and substituting the geostrophic divergence for the divergence in equa- 
tion (2.1), the vorticity equation simplifies to 



%€&) = ° 



which expanded becomes 



DJl 

Dt 



£ 




(2.4) 



Equation (2.4) shows that parcels moving poleward will increase 
their relative vorticity and vice versa. Expanding (2.4) further gives 

( 2 . 5 ) 



where 



3^ = 1 



( 2 . 6 ) 
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Substituting (2.6) into (2.5) gives the following 

' 7 ft 

(A) (8) LO (.O') 

Equation (2.7) is a prognostic equation with a single variable, z. 
Term A is the local change, term B has the mathematical form of the 
Helmholtz term but with variable coefficients, term C is the advective 
term, and term D is the ^ term whose sign and magnitude now depend on 
the relative vorticity. 

It is proposed to test equation (2.7) with and without term B and 
with a tuning constant on term D. 

A basic linearized expression can be derived for equation (2.3). 
Assuming only a basic zonal current U and that products of perturba- 
tions are small and can be ignored, equation (2.3) becomes 

H + U t* + V <3 = -U-' P . (2.8) 



Now assuming a solution of the form 

» a (JL - Ct/ ) 

\ 

equation (2.8) reduces to 



e 



C/ - XT . 

This shows that waves will travel with the speed of the basic 
current. This analysis shows that this model should move waves faster 
than models whose waves move with the speed of the Rossby waves [3], 

6 ■ U - • 
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3. 



Procedure. 



Data for 850, 300, and 500 mb were obtained from Fleet Numerical 
Weather Facility, Monterey, California, for summer and winter situations. 
These, data had been smoothed using a 15- pass, 5-point smoother over a 
3969-point grid (grid distance equal to 381 km at 60N) using equation (3.1). 

S - D + X V 2 D , (3.D 

where k is a constant smoothing factor. 

Using this smoothed field, equation (2.7) becomes 

V 2 ~ V J~t; ' 

which after applying finite-difference approximations [5], becomes 

C = - >V <j( (3.2) 

where ™ ^ ^ S ’ ^ J/rvP 3 * 

Equation (3.2) was programmed for the CDC 1604 computer using a 
symbolic coded relocatable assembly program (SCRAP). 

The equation was solved for A S and then time stepped in increments 
of one hour using a forward time difference (3.3) the first hour and 
centered time differences (3.4) for succeeding hours. 

St+i = S-t-v AS , 0.3) 

S-t+i = St-\ \ 2 . A S , ( 3 . 4 ) 

At the end of 24 hours a prognostic map was plotted along with the 
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current analysis. From these two fields an error field was determined 
and from this error field a pillow (3.5) and RMSE (3.6)were computed for 
all points north of the equator. 



Pillow 

RMSE 



= 4 (A-B) 
N 




H 



(3.5) 

(3.6) 



A represents the prognostic field, B the verifying field, and N the 
number of points used. 



16 



4. 



Results . 



Table 4.1 lists the pillow and the RMSE for tests run with the B 
term of equation (2.7) in and out for summer and winter situations. 

All runs were for 24 hours. This table shows that the B term tends to 
increase the error. 

Figure 4.1 is the prognostic map with the B term in. Superimposed 
on this map in dashed lines are the major troughs and centers of the 
prognostic map with the B term out (fig. 4. 2). Analysis of this map 
shows that the B term does not affect the movement of the waves but 
rather the intensity of the systems. The B term deepens lows and troughs 
and builds highs and ridges. The high over Northern Russia is 5349m 
with the B term in and 5321m with the B term out. It verifies as 5250m. 
The low over Greenland, Northern Canada, and the North Pacific covers 
a much broader area and is more intense with the B term in than the same 
low produced with the B term out. The verifying low is in better agree- 
ment with the low shown on the map with the B term out. 

Another matter which must be considered is the solution time on the 
computer. With the B term in, it takes 105 seconds to solve the equa- 
tion for one hour, 90 seconds of which is used to relax the "Helmholtz" 
term while it takes only 35 seconds to solve the equation with the B 
term out, 20 seconds of which is used in relaxation. Time on a computer 
is a premium and must be considered in any model. 

Figure 4.2 is the prognostic map with the B term out. Superimposed 
on this map are the major troughs and centers of the verifying map 
(fig.4.3)» Analysis shows that this model moves the waves faster than 
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B in vs B out 
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24-hour forecast for 500 mb 



from 8 December 1965, 



Figure 



4.1 
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B term in 




20 



their verifying speeds. Centers also move too fast and are in general 
too intense. 

Recalling that in the perturbation equation (2.8), the ^ terms 
cancelled each other and thus resulted in the waves moving with a speed 
equal to the zonal wind, it was felt that by changing the sign of the ^ 
term in our model (D term) we would put some of the ^3 effect back in 
and slow the waves down. 

Table 4.2 lists the error fields for summer and winter situations 
at 500 mb. This table shows that with a tuning constant of minus 4 on 
the D term, the model produces the best results. 

Figure 4.4 is the prognostic map with the B term out and a tuning 
constant of minus 4 on the D term. Superimposed in dashed lines are 
the major troughs and centers of the verifying map. 

Analysis of this map shows that th.e waves have slowed and now move 
with a speed consistent with the verifying waves. The centers appear 
to be in quite good agreement both in location and in intensity. It was 
found that increasing the tuning constant resulted in continued slowing 
of the waves. 

It must be remembered that the verifying maps are prejudiced by 
FNWF's prognostic model since this is used as a first guess to the verify 
ing data. Therefore in sparse data areas (mostly open oceans) and wher- 
ever questionable data are received, the prognostic data are used. 

Figure 4.5 is a 500-mb prognostic map for summer with the B terra 
out and a tuning constant of minus 4 on the D term. Figure 4.6 is the 
verifying map. Of interest in comparing these maps is the agreement on 
the trough coming out of the 5426m low east of Greenland and extending 
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B term out 
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Figure 



4.4 



24-hour forecast for 500 mb from 8 December 1965, B term out; 
tuning constant minus 4 
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Figure 



4*5 



24-hour forecast for 500 mb from 3 
tuning constant minus 4 



July 1966, B term out; 
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Figure 4.6 500-mb analysis, 4 July 1966 
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down the European coast and out into the Atlantic west of Gibraltar. 
Again a trough in Western Russia which extends over the Black Sea and 
into the Eastern Mediterranean compares favorably. The low over Japan 
agrees in position but is 20m higher than the verifying height. The 
trough coming out of Alaska and down Western Canada into the United 
States does not show a closed low over North America as does the pro- 
gnostic map but a low of 5659m is indicated on the verifying map and 
could possibly have been drawn as a closed system. 

Table 4.3 lists the error fields for the 850-mb level for both 
summer and winter situations. This table shows a substantial decrease 
in the error values. Since this model is divergent, it might be expect- 
ed to perform better at a level where divergence has more effect. It 
should also be remembered that the gradients are not as strong as at 
500mb. 

At 850mb there does not seem to be any preference for a particular 
tuning constant. Figures 4.7, 4.8, and 4.9 are 24-hour forecasts made 
from 10 July 1966 with tuning constants of minus 1, minus 2, and minus 
4, respectively. Investigation of these maps will show no significant 
difference in location of the troughs and ridges or in intensity of the 
systems. Table 4.3 shows that varying the tuning constant does not ap- 
preciably change the RMSE value for any given day. 

Figure 4.10 shows a typical prognostic map for 850mb and figure 
4.11 is the verifying map. Investigation of these maps will show quite 
good agreement in both the location and intensity of the major systems. 

Test results, error fields, and investigations of the other maps 
in both summer and winter show similar patterns and substantiate these 
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Figure 4.7 24-hour forecast for 850 mb from 10 July 1966, B term out; 
tuning constant minus 1 
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Figure 4.8 24-hour forecast for 850 mb from 10 July 1966, B term out; 
tuning constant minus 2 
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Figure 4.9 24-hour forecast for 850 mb from 10 July 1966, B term out; 
tuning constant minus 4 
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Figure 4.10 



24-hour forecast for 850 mb from 8 December 1965, B term out 
tuning constant minus 4 
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Figure 4.11 850-mb analysis, 9 December 1965 
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initial results 



Table 4.5 is a comparison between the error fields obtained with 
this model and the error fields of FNWF*s model for the same dates at 
500 mb. The table shows that this model produces results comparable to 
those of the present operational model and in most cases might surpass 
it if this model were used as a first guess field for analysis. 

Since FNWF did not have an error field for 850mb, it was necessary 
to look elsewhere for comparison purposes. In 1962, Nicholson [4] in- 
vestigated a divergent model proposed by 'Arnason [1]. On the basis of 
the good results obtained from initial testing, ‘Arnason* s model was 
recommended for further investigation. Although an exact comparison is 
not possible since the days chosen for testing were not the same, it is 
felt that some benefit can be derived by comparing the order of magnitude 
of the error fields for each model. 

Table 4.6 3hows the 850-mb RMSE values for both summer and winter 
situations. The verification scores for ‘Arnason* s model are those 
obtained using the optimum version of the model. The values entered for 
Duthie*s model are those showing the largest errors. 

During initial testing at 300 mb, it was found that a relaxation 
coefficient of 17cm was to small. This coefficient required the computer 
to make in excess of fifty passes to arrive at a solution. Since each 
pass takes approximately a second and a half, it required close to one 
hour to make a 24-hour forecast including printing. It was decided to 
try 35cm on a trial basis. This cut the running time in half and was 
used throughout the remaining tests. However it is felt that this 
coefficient could still be increased and further experimentation is 
recommended. 
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Table 4.4 lists the RMSE values for tests run at 300 mb for both 



summer and winter situations. As was expected, it shows that the 
errors have increased substantially, being more than double those for 
500 mb. This can partially be explained by the large gradients found 
at this level. 

Figure 4.12 is the 24-hour prognostic map from 11 December 1965 with 
the B term out. Figure 4.13 is the verifying map. Investigation re- 
veals that the model is moving waves too rapidly eastward. Figure 4.14 
is the 24-hour prognostic map for the same time but with a tuning con- 
stant of minus one on the D term. Comparing figure 4.14 with figure 
4.12 shows that the tuning constant has slowed the waves down somewhat 
but they are still not in agreement with their verifying positions. 
Analysis of other maps with various tuning constants shows that none 
was able to compensate for the strong zonal wind at 300 mb. Figure 
4.15 is the prognostic map for a summer situation with the B term out 
and a tuning constant of minus one on the D term. Figure 4.16 is the 
verifying map. In summer when the zonal wind is not quite as strong, 
the errors are reduced but the forecast waves are still moving too 
rapidly. 

It is felt that by using the gradient wind instead of the geostrophic 
wind in the C term of the model, it will reduce the speed of movement. 
Since the patterns of the prognostic and verifying maps are quite 
similar this modification should reduce the overall error field. 

During the course of this investigation it was found that the use of 
RMSE to define an error field did not do an adequate job. RMSE is only 
a measure of the difference in numerical values of two parameters. It 
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Model vs Fleet Numerical 
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Figure 4.12 24-hour forecast for 300 mb from 11 December 1965, 
d term out 
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Figure 4.13 300-mb analysis, 12 December 1965 
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Figure 4.14 24-hour forecast for 300 mb from 11 December 1965, 
B term out; tuning constant minus 1 
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Figure 4.15 24-hour forecast for 300 mb from 3 July 1966, 
B term out; tuning constant minus 1 
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Figure 4.16 



300-mb analysis, 4 July 1966 
43 



gives no indication as to how the maps resemble each other. It is felt 
that using RMSE in conjunction with a correlation coefficient, which 
is a measure of how much the products look alike, would provide a more 
useful error field. For example, a 300-mb map compared to a 500-mb map 
would show a high correlation coefficient but would also have a high 
RMSE. A high correlation coefficient together with a low RMSE would 
indicate that the prognostic map not only looked like the verifying 
map but also had small numerical errors. 
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5. 



Conclusions . 



On the basis of the results of this investigation, the model 
produces a prognostic map which is a very good representation of the 
verifying data. At both 500 mb and 850 mb the waves of the 'tuned' 
model move with a speed consistent with the verifying waves and the 
location and intensity of the centers agree remarkably well. At 300 
mb the model tends to move systems too far. This is due to the use of 
a geostrophic wind with the strong zonal components at this level. 

In view of the close agreement at 500 mb and below, it is felt 
that further investigation is warranted, incorporating the gradient 
wind in lieu of the geostrophic wind, the geostrophic relative vorticity 
(equation 2.6 with the second term in the brackets omitted) and daily 
testing begun over an extended period to determine the feasibility of 
this model for possible operational use. 
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